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Abstract: We study bound-state effects on tlie pair production of gluinos at liadron 
colliders, in a context of the minimal supersymmetric extension of the standard model. 
Due to the expected large mass and the octet color-charge of gluinos, the bound-state 
effects can be substantial at the LHC. We find significant deformation of the invariant-mass 
distributions of a gluino-pair near the mass threshold, as well as an additional correction 
to the total cross-section. Both the invariant-mass distribution and the correction to the 
total cross section depend crucially on the decay width of the gluino. 
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1. Introduction 

At the CERN Large Hadron Cohider (LHC), discovery of new heavy particles at TeV 
scale is expected as an evidence of physics beyond the Standard Model (SM). One of the 
well-investigated models beyond the SM is the minimal supersymmetric extension of the 
Standard Model (MSSM). At hadron colliders, colored particles can be produced copiously. 
In the MSSM, the colored particles are superpartners of quarks and gluons, which are called 
squarks and gluinos, respectively. The LHC has high discovery potential for those parti- 
cles, where masses up to 2 TeV can be explored [1,2]. Current limits on their mass and 
production cross-sections are set by measurements at the Fermilab Tevatron collider [3, 4] 
(See also Ref. [5]). 

In the MSSM, due to the i2-parity conservation, superpartner particles or sparticles, 
are produced in pair. The total and differential cross-sections at hadron colliders have 
been known for long time in the leading-order (LO) [6-8], as well as in the next-to-leading 
oder (NLO) including SUSY-QCD corrections [9-11]. The NLO corrections are found to 
positive and large. The origin of the large corrections comes from the threshold logs due 
to the emission of the soft /collinear gluons, and also a Coulomb singularity due to the ex- 
change of Coulomb gluons between the final-state particles. Both terms become significant 
at the partonic threshold region, where they should be considered to all-order. Recently, 
the all-order summation of the threshold logarithmic corrections has been performed for 
the sparticle-pair production processes at hadron colliders in Refs. [12-14]. 
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The aim of this paper is to examine the effects of Coulomb corrections to all-orders 
in the gluino-pair production process. Similar studies for gluino-pair and also squark-pair 
productions can be found in Ref . [14] , where the Coulomb corrections are taken into account 
to all-orders by using the Sommerfeld factor [15, 16]. However, to our knowledge, there has 
been no study including the bound-states effect which have been found significant for the 
top-quark pair production at the LHC [15, 17, 18]. ^ The bound-states of a pair of gluinos 
are called gluinoniuni [20-28]. When the Coulomb force is attractive, it causes not only the 
corrections to the cross-sections above the mass threshold [14], but also formation of the 
gluinoniums below the threshold. The bound-state formations can be taken into account 
by summing the Coulomb terms to all-orders or non-perturbatively. The Green's function 
formalism [29, 30] has been developed for this purpose, which incorporates the finite-width 
effects of the constituent particles. 

The paper follows the recent studies for the bound-state effects in the top-quark pair 
production at hadron colliders [17,18]. The invariant-mass distribution of the top-quark 
pair is distorted significantly near the threshold region [17, 18], and we may expect similar 
effects for the gluino-pair. 

The paper is organized as follows: In Sec. || and Sec. |^, we briefly review the basic prop- 
erties of the gluino decay-width and the gluino-pair production cross-section, respectively. 
Then in Sec. ^, we discuss binding corrections to the gluino-pair production. In Sec. |5|, 
we discuss effects of the initial-state radiation and the hard-vertex correction. Finally, we 
summarize our findings in Sec. ^ 



2. Gluino decay- width 

In this section and the next section, we briefly review the basic properties of the gluino 
that are relevant to our studying; the decay-width and the pair-production cross-section 
at hadron colliders, respectively. For more details, see e.g. Ref. [31]. 



The effects of gluinonium formation at hadron colliders depend strongly on the gluino 
decay width. In particular, the following regions should have qualitatively different signals: 



A 
B 
C 
D 



> \Eb\ 

\Eb\ > > Tgg 

\Eb\ > Tg > Tgg 

< ^aa 



(2.1a) 
(2.1b) 
(2.1c) 
(2.1d) 



Here, Eb < is the binding-energy of gluinoniums, and F^^ denotes the partial decay- 
width of the gluinonium annihilation into gluons. In the region A, the produced gluinos 
decay before the gluinonium formation, and the binding effects cannot be observed. In 
the region B, the binding effects are expected to enhance the pair production cross-section 



""^We are informed that the bound-state effect on the squark pair production is examined in Ref. [19]. 
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below the threshold, especially at around the location of the ground state energy, as is ex- 
pected for the top-quarks [15, 17, 18]. In the region C, a few narrow gluinonium resonances 
can be produced, whose decay may still be dominated by the constituent gluino decays. In 
the region D, the produced gluinonium will decay mainly into gluon jets, and hence will 
disappear without leaving a detectable trace in hadron collider environments. 

In Fig. ^, we plot the decay width of the 
gluino as a function of the gluino mass. The 
decay- width is calculated for g — > qq^'^W 
and g qqB decays [32], where winos {W) 
and a bino {B) are superpartners of SU(2)l 
and U(1)y gauge bosons, respectively. Solid, 
dashed and dot-dashed curves are for the 
squark mass ruq = 500 GeV, 1 TeV and 1.5 
TeV, respectively. For simplicity, we assume 
that all the super-partners of 5 light-quarks 
with both chiralities have a common mass 
rriq, and the gaugino masses satisfy the re- 
lation rrig : : = 7 : 2 : 1, which 
are valid in several SUSY breaking scenar- 
ios. Neither the contribution from the top 
squarks nor the gaugino-higgsino mixing are 
considered for brevity. 

The upper dotted-line shows the mag- 
nitude of the binding energy of the ground 
^So(l) gluinonium, which is the most deeply- 
bounded color-singlet gluinonium, see sec- 
tion ^ in detail. It is estimated by using a 
Coulombic potential V{r) = —CACts/f with the color-factor Ca = 3 for the color-singlet 
gluinoniums. Within this approximation, the binding-energy is calculated as 

Eb['So{1)] = -^mgal{f,B)MS^ (2.2) 

where we take the MS renormalization scale to satisfy fiB = CA'l^^gC^s{^J'B)M^/'^^ which 
corresponds to a half of the inverse Bohr radius of the ^S'o(l) gluinonium. It is known that 
this scale choice makes the QCD higher-order corrections to the binding energy small. For 
the gluino-mass from 200 GeV to 2 TeV, the binding energy grows from ~ 10 GeV to ~ 50 
GeV. 

The lower dotted-line shows the partial decay- width of the ground ^5'o(l) gluinonium 
annihilation into two gluons; 

r,,[^5o(l)] = 18vra^(^)5^.^, (2.3) 



500 GeV 




1500 2000 



Figure 1: The gluino decay- width as a function 
of the gluino mass, calculated from g —f qq^'^W 
and g qqB decay modes. We consider de- 
cays into 5-flavor massless quarks only where 
the corresponding squarks have a common mass 
rUq = 500 GeV, 1 TeV or 1.5 TeV, and the gaug- 
ino masses satisfy rrig : rriy^ : = 7 : 2 : 1. 
The dotted-lines show the magnitude of the 
binding energy of the ground ^S'o(l) gluinonium 
{—Eb\^Sq]), and the partial decay-width of the 
ground ^5*0(1) gluinonium annihilation into glu- 
ons [Tggi^So]). 
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Figure 2: Feynman diagrams for gg gg process and qq gg process at the tree-level. Dashed- 
lines represent squarks, (Jr or q^. 

with |V'(0)|2 = ClmlaUfiB) 

j^/Svr which is also obtained by using the Coulombic potential 
approximation. For the same gluino-mass region as above, it grows from several hundred 
MeM to one GeV. 



The three curves for Tg show steep gradient at rUg ~ rrig. Above the squark mass 
threshold, the gluinos decay into a squark and a quark and the width is proportional to 
aging. Below the two-body decay threshold, the gluino decay width drops quickly down 
to the three-body decay width with the electroweak coupling. The gluino width curves 
cross the binding-energy curve at around nig ~ l.Sniq, while they cross the gluinonium 
annihilation width just above nig = rriq. Therefore, if the gluino is lighter than all the 
squarks, the region D ( ^.ld| ) would likely be realized, and binding effects cannot be observed 
at hadron colliders. On the other hand for rrig > nig, depending on the mass ratio nig/niq, 
regions ^ to C (2.1a-c) would be realized. 



3. Gluino-pair production cross-section 

Now, we review the gluino-pair production at hadron colliders, focusing our attention on 
the threshold behavior and the color structure of the scattering amplitudes. 

At the parton level, there are two leading subprocesses for the gluino-pair production; 

g{pi,\i,ai) + g{p2,X2,a2) gipsAs^a^) + 5(^4, A4, 04), (3.1a) 
q{pi,Xi,ii) + q{p2,X2,i2) ^ giPsAs^a^) + 5(^4, A4, 04). (3.1b) 

The tree-level Feynman-diagrams are shown in Fig. ^. Here, pk and Xk are momenta and 
helicities of particles, and are color indices of gluons (gluinos) and quarks, respectively. 
We normalize Xk to take ±1 both for fermions and gluons. At the LHC, due to the steep rise 
of the gluon density in the small momentum- fraction region, gluon- fusion process ( 3.1a| ) 



dominates the cross section for nig < 1 TeV [9]. For nig > 1 TeV, the qq annihilation 
process ( p.lbj ) can contribute significantly to the total cross-section. 

3.1 gg gg process 

First, we review the gluon- fusion process. We present helicity amplitudes for this process 
by decomposing the color factors into symmetric and antisymmetric parts: 

Mgg-.~gg{p^; Xf, a,) = [1 {F'^^ F^^} ^^^^ Mggipf, A,) + ^ [F'^^F-%^^^ Ngg{p,; A,)] • (3.2) 
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Mggipk;X,-X,\,~X) = -~X (3.3c) 
M,,(p,; A, -A, A, -A) = -2 'T^'Z17r' , (3.3d) 



Here gs = v'SttoJ is the QCD couphng constant. = —if^^^ is the color matrix in the 
adjoint representation. The amplitudes can be obtained by replacing only the color-factors 
from those of the heavy-quark pair production, gg QQ. The color-symmetric part and 
antisymmetric part of the amplitudes are summarized as follows; 

Mggipk] A, A, A, A) = ^J^^-^_^^tco^2g ' (3.3a) 

Mgg{pk;X,X,~X,-X) = 0, (3.3b) 

- Arrig [5 sir? 6 

1-/3^ COS^ I 

^/3sin6'(AA + cos6') 
1 - /?2 cos2 e 

for the symmetric part, and 

Nggipk-Ak) = cose ■ Mgg{pk;Xk). (3.4) 

for the anti-symmetric part. Here, is the angle between the 3-momenta pi and ps, and 
13 = ^Jl — im'^/s is the velocity of the gluino in the partonic center-of-mass frame with 

S = {pi +P2)^. 

The product of the two color-octet states transforms as 808 = 1© 8s© 8a © 10 © 
gg amplitudes, the color basis for these states can be taken as [12]; 

ijaia2^a3a4^ (3.5a) 
3^aia2fe^a3a46^ (3.5b) 

5 

Ijaia2feja3a4fe^ (3.5c) 
1 1 

_ ^^01113^0204 _ ^0104^0203 -J ^01026^0304?) ^-g g^-j 

1 13 

_ ^^0-113^0204 _|_ ^0104^0203 -J _ _^Oi02^0304 _ _^Oi026^03046 ^2 Qg-j 

2 8 5' 

where the bases ci to C5 are for 1, 8s, 8a, 10 © 10, 27, respectively. The color factors in 
Eq. (|3.2| ) are then expressed as 

1 {^01 ^02 1^^^^ ^ 1 1^^01036^02046 ^ ^01046^02036^ ^ ^ 3^^ _ 



27. 


In the 


gg 




1 : 


Cl 




8s : 


C2 




8a : 


C3 


10 ( 


BIO : 


C4 




27 : 


C5 



^ [^"^^"']o304 = I [n^'T'''' - r^-^V"^"^') = \c^- m 

The C4 term for 10 © 10 state does not appear at the tree-level. Near the threshold, only 
the color-symmetric part of the amplitude Mgg survives, and the anti-symmetric part Ngg 
is suppressed by j3. Thus, production cross-section for the color-8A state is suppressed 
by 0{(3'^) near the threshold. All the other color-states have the same production ampli- 
tude at the tree-level, with the normalization constant (|3cip, |3/2c2p, |c5p) = (9, 18,27). 
Therefore, 1/6 of gluino-pair produced near the threshold is in the color-singlet, while 1/3, 
1/2 are in the 8s, 27 color-states, respectively. 
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3.2 qq gg process 

Next, we examine the qq annihilation process. In contrast to the qq — > QQ process where 
only the color-octet state are produced, the color-singlet gg state can also be produced 
due to the squark-exchange diagrams; see Fig. 0. 



Helicity amplitudes for this process are written as; 

1 



9s 



The color-symmetric amplitudes are 



Mj9(pfc;A, -A, A,A) = A 



Arrig (3 cos 6 sin 9 



Mqq{pk;X,-\,X,-\) = -2/?(l + AAcos( 



and the anti-symmetric ones are 



^qq(Pfc;A, -A,A,A) = A 



4m n 



■ smf 



Ax — XX cos 9 
^Al-(3^ cos2 0' 



Ax 



(3.8) 

(3.9a) 
(3.9b) 



42 



NqgiPk] A, -A, A, -A) = 2 (AA + cos 9) 



- /32 cos2 9 
Ax - AA^2 (,os 
A\ - /?2 cos2 9 



(3.10a) 
(3.10b) 



Here, we define Ay 



)/s where g+ = qR and g_ = qi. Due to the 



U = 1 - 2(m2 - 

chirality conservation, amplitudes for Ai = A2 vanish and only the squark with the chirality 
of the incoming quark A = Ai contributes. In the threshold limit /3 — > 0, the color- 
symmetric amplitudes behave as Mqq ~ 0{l3) and the anti-symmetric ones as Nqq oc 
(m| — m?^)/(m| + m|^). Note that, the first term in the square bracket in Eqs. ( |3.10| ) comes 
from the s-channel gluon exchange diagram, while the other term comes from squark- 
exchange diagrams. When Ax = 1 i.e. when the gluino mass and the squark mass are 
degenerate, the s-channel diagram and squark-exchange diagrams interfere destructively [9] 
to make the full amplitudes (|3.8[) color-symmetric. 

The color states of the produced gluino pair are found from 



1 1 



1 



(3.11a) 
(3.11b) 



The first term in the right-hand side in Eq. ( 3.11a| ) represents the color-singlet, and the 
others are color-octets either symmetric ( [3.11a ) or anti-symmetric ( [j.llbl) . Thus, in qq 
annihilation, production of the color-singlet state is suppressed by 0{j3'^) and the color- 
octet state (8a) dominates near the threshold, unless the gluino mass and the squark mass 
are degenerate. 
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4. Binding corrections 



In this section, we describe the bound-state effects in the gluino-pair production. The 
partonic cross-section including bound-state correction is expressed as; 



(c), 



Im[GW(0,^ + ir^)] 
Im[Go{0,E + irg)] ' 



(4.1) 



where al q is the tree-level cross-section. Here, we specify the color-state of the gluino-pair 
denoted by c, and i represents the initial partons i = gg or qq. The Green's function G^^^ 
is obtained by solving the Schodingcr equation; 

{E + iTg)-\-— + y (r) 11 G^'^) (x; E + iTg) = S^{x), (4.2) 

where E is an energy of the gluino-pair measured from their mass threshold. We define it 



as. 



E 



M - 2mg (M < 2mg) 
m~g0^ {M>2m~g) ' 



(4.3) 



where M is the invariant-mass of the gluino-pair. V'^^\r) is QCD potentials between the 
gluino-pair depending on their colors. Go is obtained by taking V^'^^ = 0, which is simply 
the tree-level Green's function of the gluinos with its finite width. 
We use the NLO QCD potential [33]; 



F(c)(r) = g('=) 



as in 



-bJms 



1 + ^ {2/3o [In ifiBT) + je] + ai} + o({^f)] (4.4) 



where the color-factor C^'^^ is known [23] to be 



{-Ca, ~2^a, ~2^a, 0, 1}, 



(4.5) 



for c = {1, 8s, 8a, 10 © 10, 27}, respectively. Thus, QCD potential for 1, 8s and 8a color- 
states is attractive, while that for 27 color-state is repulsive. 7e = 0.5772 . . . denotes the 
Euler constant, and the other constants are; 



11 



31 



/?o = -w^A - 7;nn, ai = —Ca 



10 



-n„ 



(4.6) 



We take the number of light-quark flavor to Uq = 5, and the scale as a half of the 
inverse Bohr radius of the color-singlet gluinonium, hb = l/2rB = CAmgasinB)/^- For 
example for rrig = 608 GeV, we obtain hb — 58 GeV and as{nB) — 0.127. With this scale 
choice, the higher-order corrections to the binding energy is small. On the other hand, the 
magnitude of the Green's function is reduced by several tens percent in an entire range of 
E, when one includes the NLO correction. 
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For the attractive force, the gluino-pair forms a bound-state with a definite spin (S) 
and an orbital angular momenta (L) [23]. Within the MSSM the gluino is a Majorana 
fermion, and the wave-function of the gg state must be anti-symmetric under the exchange 
of the two gluinos. This leads to the restriction on the spin and orbital angular momenta 
of the gluinonium, so that the color-symmetric pair in (1, 8s, 27) forms only even L + S 
states, while the color-antisymmetric pair (8a) forms only odd L + S states. Possible 
bound-states for each color-state are summarized in Table |l|. We also list the non-zero 
components in the tree-level amplitude for each process. Blank states in i = gg and i = qq 
rows are missing at the tree-level. For example, due to the Yang's theorem, color-symmetric 
J = 1 states are forbidden in gluon-fusion process. Likewise, = or J = states are 
forbidden in i = qq, due to the chirality conservation. 

For the repulsive force, the gluino-pair does not form a bound-state, but the Green's 
function formula is still useful for summing over Coulombic corrections. 

In our analysis, we include only the S-wave Green's function. The other states are left 
free, since they are not affected significantly by the binding effects. The extraction of the 
certain partial-wave amplitude is straightforward from the helicity amplitudes given in the 
previous section. 

The gg invariant-mass distribution in the hadronic collisions is given as [17, 18]; 

where r = M^/s. The partonic luminosity is defined by 

^(r;^F)= J dxi J dx2fa{xi, ^F)fbix2, IJ'f) 6 {t - X1X2), (4.8) 

{a,b} 

where the summation is over {a,b} = {g,g} for i = gg, and {a,b} = {q,q},{q,q} with 
q = u,d,s,c,b for i = qq. We use 5-flavor parton distribution functions, and the 5-flavor 
strong coupling constant even for > mt. 

The convolution formula in Eq. ( [4. 71 ) is based on an assumption of the factorization of 
the soft /collinear gluon emission and the Coulomb corrections^. Thus, the soft/collinear 
gluon emission is described by the ISR (initial-state radiation) functions Fj^^\ and the 



^Recent paper by M. Beneke, P. Falgari and C. Schwirin shows the factorization property of the Coulomb- 
gluon and the soft-gluon [34]. 



color 


symmetric (l, 8s, 27) 


anti-symmetric (8a) 


gg 


■^Sq, ^-Po,i,2, 'D2, ■ ■ ■ 




i = gg 


■^Sq, ^Pq, 2, 'D2, ■ ■ ■ 


'Pi, '^1, 3, ••• 


i = qq 




'Si, 3^1,2,3, ••• 



Table 1: Summary of gluinonium states, and possible tree-level production from i = gg and i = qq. 
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Figure 3: gg invariant-mass distribution at the LHC, for rrig — 608 GeV and Tg = 5.5 GeV 
(SPSla). Result in the Born level (dotted), 0{as) Coulomb correction (dashed), and the all-order 
Coulomb correction (solid) are plotted. Right figure is the same, but enlarged for the threshold 
region. All-order Coulomb correction for individual 1 and 8s color-states in the gluon-fusion process 
are also plotted in dot-dashed lines. 



Coulomb corrections are included in the partonic cross-section d\ , as in Eq. ( [4.1[ ). 
represents the hard-vertex correction, which is a process-dependent constant. 

In this section, we neglect the ISR effect for simplicity, and set F- '^\z) = 6{1 — z) with 

(c) 

=1. The effect is discussed in the next section. 

In the left figure in Fig. ^, we show the invariant-mass distribution of the gluino-pair 
including bound-state corrections, for the LHC with ^/s = 14 TeV. The results are given 
for rrig = 608 GeV, a la SPSla [35], and for Tg = 5.5 GeV for nig = 547 GeV. We use 
the CTEQ6L PDF parameterization [36], and set scales as fiji = fip = mg- Predictions 
in Born-level, including 0{as) Coulomb correction, and all-order Coulomb correction are 
plotted in dotted, dashed and solid lines, respectively. Compared to the Born-level result, 
the 0{as) correction as well as the all-order Coulomb corrections bring large enhancement 
near the mass threshold (M ~ 2mg). Moreover, the distribution for the all-order correction 
indicates the resonance peaks below the threshold. Due to the binding correction and the 
finite- width effect, the gluino-pair invariant-mass distribution arises rapidly below the mass 
threshold. Thus the effective pair-production threshold is not at 2mg but at M ~ 2mg — 30 
GeV. The order of the deviation is about the binding-energy of the ^5*0(1) state, thus 
0{mga1). On the other hand, the enhancement disappears in the high invariant-mass 
region. 

The all-order Coulomb correction enhance the total cross-section by about 15% from 
the cross-section in Born level, where a substantial portion from below the threshold region 
contributes. We discuss this result in more detail below. 

The right figure in Fig. ^ is the same, but enlarged on the threshold region. For the 
all-order correction, individual contributions of the color 1, Sg states are also plotted. But 
those for other color-states, 8a, 27 which have no enhancement near the threshold are 
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not shown, qq annihilation process is neghgible in this case. In the ah-order distribu- 
tion, two resonance peaks can be found at ~ —22 GeV and E ~ —7 GeV. The first 
peak corresponds to the ground state of the color-singlet ^Sq resonance, and the second 
peak consists of the second excitation of the ^5*0(1) resonance and the ground state of the 
-"^ 50(88) resonance. 

Fig. ^ is also the invariant-mass distribu- 
tions at the very threshold region, considering 
a variety of the gluino decay width. Results 
for rrig = 608 GeV, for F^ = 20 GeV, 5.5 
GeV, 1.5 GeV and 0.5 GeV are plotted in the 
dotted, dot-dashed, dashed and solid lines, re- 
spectively. 

According to Fig. |, F^ = 20 GeV is the 
same order as the gluinonium binding energy 
I^bI ~ 22 GeV for rUg = 608 GeV. Thus it 
corresponds to the region A according to the 
classification defined in Sec. In this case, 
the invariant-mass distribution shows no res- 
onance peaks, but a gradual slope until rather 
below the threshold. This is due to the finite 
width effect, giving the production probabil- 
ity of the events with a off-shell gluino. The 
Coulomb corrections act to enhance the pro- 
duction ratio at the threshold region. Next, F^ = 5.5 GeV is the same value as used in 
Fig. ^ which corresponds to the region B. This case, the distribution shows a broad reso- 
nance at ~ —22 GeV, corresponding to the ground ^5o(l) gluinonium. The resonance 
for the ground ^50(88) gluinonium at £^ ~ —7 GeV is barely seen. In the distribution for 
F^ = 1.5 GeV (region C), the two separated resonances for the ground ^50(1) and ^S'o(8s) 
gluinoniums can be seen. Finally, for F^ = 0.5 GeV (region D), several resonances can be 
found in M < 2mg. Actually, the first peak for the ground ^S'o(l) resonance is quite sharp. 
The second peak at E' ~ —7 GeV which consists of the ground ^50(88) gluinonium and 
the second excited ^5*0(1) gluinonium is also clearly seen. 

As mentioned in Sec. ^, for the regions C and D, the partial decay width of the gluino- 
nium annihilation into gluons can not be negligible. For nig = 608 GeV, the two gluon 
decay-width of the ^5'o(l) gluinonium is about Fgg[^5o(l)] ^ 0.75 GeV, see Eq. ( |2.3| ) and 
Fig. ||. Note that the gluinoniums have partial decay width 2Tg, due to the decay width 
of the constituent gluinos. Thus, the branching ratio to the "hidden" gluino decay may be 
estimated by B{{gg) — > gg) = Tgg/{2Vg + Tgg). Using this estimation, about 20% (40%) of 
the resonance events for F^ = 1.5 GeV (0.5 GeV) in Fig. ^ would decay into jets. Therefore, 
only 80% (60%) of the resonance events would leave detectable signals of the gluinos at 
hadron colliders. 



0.08 - 




gluino-pair iiivarianl-mass [TeV] 

Figure 4: Invariant-mass distributions at the 
threshold region, at the LHC, ^ = 14 TeV. 
Results including the all-order Coulomb cor- 
rection, as well as no Coulomb correction are 
plotted for m-g = 608 GeV, for Tg = 20 
GeV (dotted), 5.5 GeV (dot-dashed) 1.5 GeV 
(dashed) and 0.5 GeV (solid). 
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A correction to the total cross-section is also an important consequence of the inclusion 
of the binding corrections. Above the threshold, the all-order summation of the Coulomb 
corrections gives the Sommerfeld factor. It enhances the total cross-section from the Born- 
level results by 10% to 20% for rrig = 200 GeV to 2 TeV, almost independent of the gluino 
decay-width. This result is rather smaller than that in Ref. [14] , due to the NLO correction 
to the QCD potential which reduces the magnitude of the Green's function even for £^ ^ 0. 

In addition, there comes a substantial portion to the total cross section from the 
M < 2mg region. The cross section in this region emerges by (1) the smearing effect due 
to the finite width of the gluino, and (2) the bound-state contributions. The former is 
proportional to Tg, while the later is almost independent of Tg. 

In Table ^, we calculate the proportion of the cross section in M < 2mg to the total 
cross-section. We define the total cross-section by integrating the differential cross-section 
over the invariant-mass from M = 2mg + Eb[^ So{l)] — 3Tg to ^/s, where Eb^^ So{l)] is the 
binding energy of the ^So{l) gluinonium. Tg = Eb[^So{1)], Eb[^So{1)]/2, 2Tgg[^So{l)] 
and Tgg[^ So{l)]/2, for nig = 200 GeV to 2 TeV are examined. The four patterns of Tg 
represent typical values in the regions A to D. In the square brackets are the same ratio 
but without the binding correction, thus these illustrate only the finite-width effect of the 
gluino. 

For Tg = Eb (region A), the ratio grows from 8% for smaller rrig to 11% for larger 
rrig , where a more than half of the contribution comes from the smearing effect due to the 
large decay width. The smearing effect decreases with Tg, and is negligible for the regions 
C and D. Then the bound-state effect which tends from 4% to 6% dominates the ratio. 
Basically, the ratio grows with nig. This is because, for larger nig, the gluon luminosity 
becomes more steep in M, thus the cross section in lower invariant-mass region is more 
significant. 

5. ISR corrections 

In this section, we discuss the effects of the ISR and the hard correction, which are the 
remaining major corrections in QCD. We treat the ISR function to 0{as) in a soft ap- 



m-g 


A: Tg = Eb 


B: Tg = Eb/2 


C: Tg — 2Tgg 


D- r-g = Tgg/2 


200 [GeV] 


7.5 [4.5] 


5.0 [1.8] 


4.0 [0.3] 


3.9 [0.1] 


400 [GeV] 


7.1 [4.2] 


4.8 [1.7] 


3.8 [0.2] 


3.8 [0.1] 


600 [GeV] 


7.2 [4.2] 


5.0 [1.7] 


3.9 [0.2] 


4.2 [0.0] 


1 [TeV] 


7.9 [4.6] 


5.5 [1.8] 


4.3 [0.2] 


4.4 [0.0] 


1.5 [TeV] 


9.2 [5.3] 


6.3 [2.1] 


5.0 [0.2] 


5.1 [0.0] 


2 [TeV] 


10.7 [6.3] 


7.4 [2.5] 


5.9 [0.2] 


5.9 [0.0] 



Table 2: Table of a proportion of the cross section in M < 2mg to the total cross-section (%), at 
the LHC with ^/s = 14 TeV. The same ratio but without the binding correction are also listed in 
the squark bracket. 
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proximation, as in Ref. [17]. The ISR function defined in Eq. ( [4.7D is written up to NLO 
as; 



with 



ln(l - z] 
l-z 



l-z 



and 



rrir. 



In 



2mn 



J1F_ 
2m?, 



1 - z 



(5.1) 
(5.2) 
(5.3) 



/3o „ 



The constants are given as [12]; 

= Ca, Aq = Cp, Bg 
Di = 0, Dss = -3, I^Sa = -3, 1^27 = -8. (5.4) 

The hard- vertex factor is obtained by matching with the fuh 0{as) correction. However, 
to our knowledge, 0{as) correction to the gluino-pair production has been calculated nu- 
merically and only for the color-summed cross-section in Ref. [9]. Accordingly, we cannot 
fixed the color-dependent factors, but only the color-averaged one. We parametrize the 
hard-vertex factors as 



with 



vr \mgj 



m-gj \2mg 



(5.5) 



(5.6) 



Extracting from the graphs in the figure 9 (a) and (b) in Ref. [9] or using Ref. [10], we find 



ha 



-0.26, 



hq = 1.5, 



(5.7) 
(5.8) 

where the later is obtained for niq/mg = 1.4, and we neglect the mass-ratio dependence. 

(c) (c) 

There is also an ambiguity in separating the non-logarithmic part in A;- and h\ . We take 

(c) 

sl scheme such that the non-logarithmic part in /c- ' is zero. 



Now, we examine the ISR effect to the differential cross section at the bound-state 
region. We calculate the ratio of da/dM at M = 2mg + Eb including the ISR to that 
without the ISR, K = o"isr/o"noISRi where da/dM with the ISR are evaluated using the 
CTEQ6M PDFs with ^i = mg. 

We also examine the ambiguity by the scale choice of /i. We find the maximum and 
the minimum of da/dM at M = 2mg + Eb during the range from /i = rrig to 2mp, and 
define the uncertainty for the cross sections without and with the ISR as 

<^{NoISR,ISR} = — 7 r- (ti.yj 

cr(/x = m-g) 
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In Table ^, we list the ratio K and the scale un- 
certainties 5 for the cross-section without and with 
the ISR. In the calculation, the decay width is set 
to Tg = however the width dependence is 

small. The ratio K is 1.2 for nig = 200 GeV and 
increases with nig to 2.8 for nig = 2 TeV, due to the 
terms including the plus distribution in i^^ , which 
take large values when s is close to M^. The scale 
uncertainty of the cross section without the ISR is 
21-25%, however that with the ISR is remarkably 
reduced to a few percent. 

In the invariant-mass distribution, the size of 
the ISR correction is not actually a constant, but 
moderately increase with M. As the result, the pro- 
portion of the cross section in M < 2mg to the total 
cross section shown in Table. |^ is slightly reduced. 

We note several remarks on the ISR. First, it is pointed out in the top-quark produc- 
tion [18] that the emission of the non-collinear gluon can give the same-order correction. 
This can be implemented into our analysis by utilizing the NLO calculation of the gluino- 
nium production which has not been available yet. Moreover, the color-dependence of 
the hard-vertex factor is neglected in our analysis, which is also resolved by the above 
calculation or the NLO calculation of the gluino-pair production with an explicit color 
projection. Due to the lack of them, our results have ambiguity in the overall magnitude of 
the ISR correction. However, we expect that these effects would not considerably deform 
the invariant-mass distribution. 

6. Summary 

In this paper, we study the bound-state corrections in the gluino-pair production at hadron 
colliders. We perform the all-order summation of the Coulomb corrections using the Green's 
function formalism, and find sizable corrections in the gluino-pair invariant-mass distribu- 
tion as well as in the total cross-section. 

The consequence of the bound-state correction crucially depends on the decay width 
of the gluino. When nig > niq, the gluino decay-width is typically between a few and 
a few hundred GeV; see Fig. In this case, the invariant-mass distribution below and 
near the mass threshold region receive a large enhancement, due to the formation of the 
gluinonium bound-states. We show that the gluino-pair production cross section below the 
threshold (M < 2nig) grows from 4% to 6% of the total cross section due to the bound- 
states formation, when the gluino mass increases from 200 GeV to 2 TeV. For the large 
case, although the resonance peaks are smeared-out, the ratio grows to 8% to 11% due 
to the smearing effect. When F^ is significantly smaller than the binding energy of the 



m-g 


K 


(^NoISR 




200 [GeV] 


1.17 


0.21 


0.05 


400 [GeV] 


1.36 


0.22 


0.02 


600 [GeV] 


1.52 


0.23 


0.01 


1 [TeV] 


1.82 


0.24 


0.02 


1.5 [TeV] 


2.25 


0.25 


0.03 


2 [TeV] 


2.78 


0.25 


0.05 



Table 3: Table of the ratio of the dif- 
ferential cross section da/dM at AI = 
2mg + Eb with the ISR to that without 
the ISR, as well as the scale uncertain- 
ties 6a/a. 
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gluinoniums, which grows from ~ 10 GeV for nig ~ 200 GeV to ~ 50 GeV for rrig = 2 TeV 
(see Fig. |l|), one or two resonance peaks arise in the invariant-mass distribution. Although 
the study of the observable signals of the gluino-pair production, such as the missing pT 
and various kinematical observables, is beyond the scope of this report, those events below 
the threshold should be taken into account in the determination of the gluino mass. 

On the other hand, when rrig < rriq, the gluino decay-width is quite tiny; see Fig. In 
this case, the gluinonium spectra are very sharp and more than two resonance peaks are 
expected (see Fig. ^. However, the produced gluinonium resonances would decay mainly 
into gluon jets and they may escape detections at hadron collider experiments. 
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